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This paper shows experimental evidence of photorefractive steady state self- focusing in InP:Fe 
for a wide range of intensities, at both 1.06 and 1.55//m. To explain those results, it is shown that 
despite the bi-polar nature of InPiFe where one photocarrier and one thermal carrier are to be 
considered, the long standing one photocarrier model for photorefractive solitons can be usefully 
applied. The relationship between the dark irradiance stemming out of this model and the known 
resonance intensity is then discussed. 

PACS numbers: 42.70.Nq, 42.65.Tg 



I. INTRODUCTION 



Photorefractive self-focusing and spatial solitons have 
been up to now the object of more than a decade 
of growing interest, both experimentally and theoreti- 
cally [U, U H, [1] • During this decade, solitary waves in- 
duced by the photorefractive effect have been shown to 
interact with each other showing specific behaviors upon 
collision Q like birth and fusion When a spatial 
solitary wave propagates in a photorefractive medium, 
it induces an optical waveguide that does not disappear 
instantly when the beam is turned off j leading to op- 
tically induced optical components such as junctions Q, 
couplers [loj and even o ptic ally induced soliton arrays 
and photonic lattices (Til. Il2l. [isj . 

The attraction towards photorefractive solitons and 
the waveguides they induce owes a lot to their poten- 
tialities in optically inducing photonic circuitry in the 
bulk of a material, on the one hand, and to their re- 
configuration properties, on the other hand. The former 
relies on the material property of retaining the waveguide 
once the beam is shut off, like, for instance, LiNbOa (l^ . 
The latter, on the contrary, relies on the possibility to 
rapidly screen out an applied field with a laser beam, to 
give rise to a waveguide, and to erase this waveguide as 
rapidly. Up to now and to our knowledge, the charac- 
teristic build-up and erase times are on the order of the 
second in the widely used Srj;Bai_j;Nb206 or sillen- 
ite Bii2Ti02o O, or up to the minute in LiNbOa [13, 
all of these experiments having been conducted at visi- 
ble wavelengths to which most of the cited materials are 
sensitive. Shorter response times around the millisecond 
have been obtained in semiconductor CdZnTe UM and 



around the nanosecond for visible ultra-high intensities 
in Bii2Ti02o M- 

If one wants to target optical networking as a poten- 
tial application of the dynamically reconfigurable pho- 
torefractive soliton technology, build-up and reconfigura- 
tion times below the millisecond is to be achieved ideally 
at infrared wavelengths compatible with optical fiber net- 
works. Thanks to their high carrier mobility and despite 
their low electro-optic coefficient, photorefractive semi- 
conductors such as InP:Fe might respond to this demand. 
Charge transport, in this material, has a bi-polar nature 
with photoexcited holes co-existing with thermally ex- 
cited electrons. This yields a peculiar resonance of the 
photorefractive Two Wave Mixing (TWM) gain as a func- 
tion of intensity [10, 

Owing to its bipolar transport character, the physics 
behind spatial soliton propagation in InP:Fc is more com- 
plicated than in other materials. It has already been 
the subject of theoretical analysis (23|. In this work, 
we present a series of recent investigations on photore- 
fractive self-focusing in InP:Fe at infrared wavelengths 
and we propose a model that takes into account the 
transport of both charge carriers. Unlike in previous 
work[23, m, [2^, the experimental results in our crys- 
tal samples and the corresponding theoretical prediction 
do not show any clear evidence of an influence of the 
intensity resonance of the bi-polar photorefractive effect, 
thus indicating that photorefractive self-focusing behaves 
to a large extent like in the case of the well-known single- 
carrier model [131. Hint is finally provided that the reason 
why qualitatively different behaviors are observed in dif- 
ferent samples might be a varying oxydo-reduction state. 
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II. EXPERIMENTAL EVIDENCE OF THE 
ROLE PLAYED BY DARK IRRADIANCE 

The bi-polar photorefractive resonance mechanism was 
first evidenced experimentally by Mainguet and Idrissi 
et. al. HO, [m and then defined by Picoli et. al. [12, 113 
as the mean intensity at which the Two Wave Mixing 
(TWM) gain is greatly enhanced. It corresponds to the 
intensity at which the intensity dependent hole genera- 
tion rate is equal to the temperature dependent electron 
generation rate. The resonance intensity is thus expected 
to depend on the temperature and has been evidenced so. 
It has also been shown that this TWM enhancement cor- 
responds to a 7r/2 phase shift between the index grating 
and the optical intensity grating, and to a corresponding 
enhancement of the index modulation amplitude. 

As a preliminary characterization of our samples, we 
evidenced and measured the resonance intensity in the 
samples used in our experiments. The iron content of 
the various samples we tested varied from 8 x lO^^cm"^ 
to 2 X lO^^cm"^, as measured by Secondary Ion Mass 
Spectroscopy. Similar behavior were observed in all sam- 
ples. The results presented in this paper were obtained 
with a 2 X lO^^cm"^ Fe doped sample. 

The resonance inte nsity has been measured with a 
standard TWM setup [lo] as described in for vari- 
ous temperatures, and particularly at 20°C, the temper- 
ature at which further self-focusing experimental results 
have been obtained. The measured resonance intensity 
at 20°C is 300mW/cm^ at a wavelength of 1.06/im and 
3W/cm^ at a wavelength of 1.56/im. This wavelength dis- 
parity is explained by the wavelength dependent photo- 
excitation cross sections. 

The experiments were conducted independently at 
both 1.06/im and 1.56/im wavelengths, the latter because 
it corresponds to a telecom wavelength and the former for 
InP:Fe large sensitivity to it. The 1.06/xm wavelength is 
delivered by a solid state diode pumped Nd:YAG laser 
whereas the 1.56/im beam is produced by a laser diode. 

The beam is focused down to a 30/im waist onto the 
entrance face of a temperature regulated InP:Fe sample. 
The beam propagates along the < 110 > direction and is 
linearly polarised along < 110 > with an electric field of 
lOkV/cm applied along < 001 >. The crystal dimensions 
are 5 x 5 x 12mm'^ along the directions < 001 >< 110 > 
and < 110 >, respectively. 

The steady state beam profile on the sample output 
face is observed through a microscope objective and a 
silicon CCD camera. The camera is sensitive enough for 
the 1.06/tm wavelength but requires to be equipped with 
a luminescent phosphorus converter to achieve 1.56/im 
observations. The observed image is then gamma cor- 
rected to account for the conversion non-linearity. 

The output beam profile is observed at steady state for 
various intensities, at a temperature of 20°C. The best 
result giving the strongest focusing is shown in figure 
[IJb) for both 1.06/(m and 1.56/tm wavelengths. Vary- 
ing the input intensity reveals that the steady-state self- 




(a) (b) (c) 

FIG. 1; (Color online) Steady state output beam profile for 
a peak intensity of 190mW/cm^ at 1.06/im (top) and a peak 
intensity of 8.7W/cm^ at 1.56/im (bottom), for no applied 
field (a) and a lOkV/cm horizontally applied field (b). The 
images total width is 200/im. The non Gaussian shape of the 
1.56/tm beam in absence of the applied field is due to the 
multimode character of the laser itself. The curves (c) are the 
horizontal beam profiles (a.u.) at beam maximum. 
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FIG. 2: Output beam width as a function of input beam peak 
intensity at 1.06/im (top) and 1.56/im (bottom) wavelengths. 
The beam diameter is normalized to the linearly diffracted 
output beam diameter. The horizontal dotted line represents 
an output beam as wide as the input one, suggesting a spatial 
soliton like propagation. 
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focusing phenomenon vanishes at low and high intensities 
and reaches an optimum focusing for an intermediate in- 
tensity, when the output beam width reaches a minimum. 

This is quantified on figure [2] where the normahzed 
horizontal output beam diameter is plotted as a function 
of the input beam peak intensity for an applied field of 
lOkV/cm. Linear diffraction is observed at steady state 
for small and large intensities. A minimum output beam 
diameter is reached for an intermediate intensity of ap- 
proximately 210mW/cm^ at l.OG^m and 7.5W/cm^ at 
1.56^m. A noticeable characteristic of both curves in fig- 
ure [2] is that the upside down bell-shape is attracted by 
the dotted line when it approaches it, i.e. when the out- 
put beam diameter approaches the input one. This evi- 
dences a domain of stability of the beam width (from 1 to 
4W/cm^ at 1.06/zm and from 5 to 50W/cm^ at 1.56/xm). 
The horizontal beam width associated to these domains 
roughly corresponds to the input beam width. This may 
suggest that the beam evolution is attracted towards a 
solitary wave solution 

We note that the width stability intensity range is 
much wider at 1.56yum than at 1.06^m: this is due to 
the fact that in the former case, the intensity that full- 
fills soliton propagation conditions is much nearer to the 
minimum of the reverse bell-shaped curve. Furthermore, 
one has to keep in mind that the beam intensity de- 
creases with propagation, owing to absorption, — 2cm~^ 
at 1.06/xm and 0.5cm~^ at 1.56^m — so that the soliton 
propagation conditions cannot be fuUfilled precisely all 
the way through the crystal, thus strengthening the need 
for stability if a spatial soliton is to propagate. 

As in many other materials, self-focusing in InP:Fc 
comes together with beam bending, as is evidenced by 
figure [T] and particularly [TJc). The steady-state beam 
shift at the output face with respect to the propagation 
axis, along the same horizontal direction as the applied 
field, is reported on figurc[31 We observe that beam bend- 
ing is negligible for sufficiently low and high intensities, 
that is, when self-focusing is absent. Furthermore, the 
maximum beam shift is obtained precisely at the inten- 
sity at which self- focusing is most efficient. 

The curves reported on figure [2] resemble the photore- 
fractive soliton existence curve — giving the soliton width 
as a function of the beam peak intensity to dark irradi- 
ance ratio — as predicted earlier [23, [23| in the frame- 
work of a single-photo-carrier model. The dark irradi- 
ance plays a central role in this theory since it is the ratio 
of the beam peak intensity to the dark irradiance which 
determines the steady state soliton width. The dark irra- 
diance is the equivalent uniform optical irradiance which 
should be shone on the crystal to generate as much pho- 
tocarriers as thermal excitation does. In other words, 
shining the dark irradiance on a given sample doubles its 
conductivity. 

If we assume that this one-photo-carrier model — 
which will be named model I in the following — describes 
our experiments in InPiFe, we can use its predictions to 
deduce the dark irradiance in InPiFe at both the wave- 
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FIG. 3: Output beam liorizontal shift as a function of input 
beam peak intensity at 1.06/im (top) and 1.56/im (bottom) 
wavelengtlis. 



lengths we have used. According to it, the minimum 
width soliton which corresponds to the maximum self- 
focusing effect is obtained for a beam peak intensity e qual 
to 3 times the dark irradiance, as confirmed recently [30| 
in a (l-l-l)D model. Using this hypothesis, the dark irra- 
diance in our sample can then be evaluated, from figure[2l 
to approximately 70mW/cm^ for the 1.06^m wavelength 
and 2.5W/cm^ for 1.56^m, at a temperature of 20°C; the 
accuracy of the value measured at 1.06/^m is low because 
of InP:Fe strong absorption at this wavelength. Let us 
note here that, in both cases, the estimated dark irradi- 
ance and the resonance intensity are on the same order 
of magnitude. 

At this stage, one may wonder why model I should de- 
scribe photorefractivity in InP:Fc since both holes and 
electrons are known to play a significant role [2^ . [2^ in 
this semiconductor — holes are photo-excited and elec- 
trons are thermally generated — whereas model I assumes 
only one carrier, both photo and thermally excited. One 
may hastily think that electrons can simply be neglected 
in our experiments. This is however not the case since 
both the dark irradiance and the resonance intensity are 
included in our measurement range: they respectively 
correspond to the equality of the concentrations on the 
one hand and of the generation rates on the the other 
hand, of holes and electrons. 

Therefore, since model I assumptions do not match 
with our experiments, we now need to ascertain by an 
independent measurement that the values of the dark 
irradiance we have deduced from figures [5] and [3] arc in- 
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deed correct. We have thus shone a uniform intensity on 
our InP:Fe sample and measured its photo-conductivity. 
This allowed us to deduce a 50mW/cm^ dark irradiance 
at 1.06/im and a 2W/cm^ one at 1.56/im, confirming the 
previous evaluation at this wavelength. The slight dis- 
crepancy is accounted for by the relative precision with 
which figure [2] allows the determination of the dark ir- 
radiance on the one hand and on the other hand by the 
difficulty to shine a strong and uniform light on the whole 
sample for the independent measurement. 

These recent experimental results evidence the fact 
that the dark irradiance plays a crucial role in pho- 
torefractive self- focusing in InP:Fe, as it does in pre- 
viously developed models, despite the greater complex- 
ity of the photorefractive mechanisms in InP:Fe. How- 
ever, as previously explained, model I does not allow 
to ex plai n our observations, nor do previously published 
ones j24 [25l. [26| . as will be detailed later on. Therefore, 
we will now tackle the theoretical problem of showing 
how and why the dark irradiance plays a significant role 
in InP:Fe. 



III. SPACE CHARGE FIELD AND DARK 
IRRADIANCE 



Let us start with the general two-carrier band 
port model [13, ll^l reduced to one dimension x, 
we will consider at steady state in the following, 
described by the set of equations below: 
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where E is the electric field, n and p are the electron and 
hole densities in the respective conduction and valence 
bands, tit = [Fe^+] is the density of ionized occupied 
traps, pt = [Fe'^+] is the density of neutral unoccupied 
traps, jn and jp are the electron and hole current den- 
sities, respectively. Nt, Njj and A^^ are the densities 
of iron atoms, the shallow donors and the shallow ac- 
ceptors, respectively. The charge mobilities are given 
by p,n for electrons and fip for holes, the electron and 
hole recombination rate are respectively c„ and Cp , T is 
the temperature and ki, is the Boltzmann constant. The 
dielectric permittivity is given by e while e is the ele- 
mentary charge. The electron and hole generation rates 



e„ and Cp depend on both thermal and optical emission 
as described by 



th 



th 



e„ + 



(2a) 
(2b) 



where the thermal contribution to the emission rate coef- 
ficient is e*'* and the optical cross section of the carriers is 
given by ct, / is the spatially dependent intensity of light. 
The approximations in both above equations are valid in 
the case of InP:Fe as stated above, with thermally excited 
electrons and photo-excited holes. The need for a large 
applied electric field leads us to neglect carrier diffusion 
when compared to drift, diffusion which is accounted for 
by all the terms which contain fc^T. 

Solving the differential system ([T]) requires defining the 
boundary conditions. To that aim, we will consider that 
the input beam is small with respect to the sample di- 
mensions and thus that its influence is limited to its sur- 
roundings. Therefore, the boundary values for £', n, p, 
"^T, Pt, jn and jp are the values reached at steady state 
in the dark. As was previously done [2O, H^, [28t . 
they can be deduced from system ([T]), setting to zero 
all derivatives and optical intensity. Following, if the 
subscript identifies the boundary condition: Eq = V/d 
where V is the applied voltage and d the crystal width; 
no = e^J^UTa/cnPTo and po = e^pPTolcpUTo where nro = 
No — Na and pTo = Nt — uto are the ionised and non- 
ionised iron densities in the dark [13, H^; jn = e//„noi5o 
and jp = ep,ppEo . 

The time evolution of the electric field E is determined 
as follows. The time derivative of Eq. (fla| is simplified 
using Eqs. (jldllellf]) to yield 



d^E 
dtdx 



-^ d{jn+jp) 

s dx 



(3) 



Making use of Eqs. (|lb|lcp and integrating with respect 
to the space variable x, it becomes 



dE 
'dt 



[(Mn" + ^J■pP) E. - (pn-na + MpPo) Eq] (4) 



At steady state, the evolution equations (jldp . ([Te]), flj) 
and (HI reduce to 



(^„n + fipp) E - {p-nTio + fipPo) Eo = 0, (5a) 

^dn 

dx 
idp 



CnUT - c„n {Nt - ut) + ^^nn^ + ^^nE^ = 0, (5b) 



dE 

ep {Nt - Ut) - CppriT - f^pP-g^ ~ ^^^'^ ^' ''^^^ 

To further simplify this differential system, the light 
intensity spatial variations will be assumed small enough 
so that the resulting free carriers spatial variations re- 
main small, and the derivatives dn/dx and dp/dx can 
be neglected in equations (|5bp and ([Sc)) . In the same 
spirit, we also assume that the field derivative dE/dx in 
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the same equations can be neglected, precisely dE/dx <^ 



The left-hand-side of Eq. (fla|) becomes zero within 
this approximation. Light intensity and temperature 
will again be considered not so large so that free carriers 
density remain small with respect to ionised iron density: 
n <^ TiT and p ^ n^, so that equation pa[) reduces to 
riT = f^TOJ and consequently from Eqs. (|5bp and (|5cp we 
get n = no and p = (o'pIpTo) I (cpriTo)- The assumptions 
made are valid for intensities small against Cpu'^Q/apPTo, 
itself greater than the TW/cm^, and at temperatures 
orders of magnitude above room temperature. Within 
these approximations, we obtain the electric field E de- 
pendency upon the local intensity /: 



E 



Irl 



Id+I 



Eo- 



(6) 



This expression is strictly the same as the one found 
before in the framework of model I [13, when diffusion 
is neglected; only the expression of the dark irradiance Id 
differs (the approximation below being valid in the case 
of InP:Fe): 



(7) 



Let us point out here that this has implications beyond 
the above equation: all the conclusions that have been 
driven, in model I, from an equation identical to ([6]) , still 
hold, if the above value of Id is accounted for. There 
is for instance no need to derive again the soliton exis- 
tence curve since the involved equations are the same. 
Therefore and using the predictions model I, our obser- 
vations of self-focusing in InP:Fe are partially explained. 
Indeed the presence of self-focusing over a broad range 
of intensity is predicted and no electron-hole resonance 
effect is observed. The beam bending is however not yet 
accounted for, which is one reason to further extend and 
refine our model. 



IV. RESONANCE INTENSITY 

As was said before, we felt the need to establish the 
theoretical analysis presented in the previous section — 
which we will now call model II — because we could not 
satisfactorily explain our observations using the theory 
available in the literature [13, [1E| for one particular rea- 
son: this theory — which we will call model III — predicts 
the existence of a resonance intensity around which the 
self-focusing phenomenon reverses itself; while we have 
successfully measured the resonance intensity in our sam- 
ples, we were never able to measure the predicted reversal 
from focusing to defocusing. 

The resonance intensity around which a self-focusing 
behavior reversal is predicted by model III is precisely 
the same resonance intensity at which TWM gain is max- 
imum. As was said before, the resonance intensity is also 



the intensity at which the generation rates of holes and 
electrons are equal: 



Ir = 



PTo 



<^pPTo 



(8) 



At resonance, the imaginary part of the space charge 
field reaches a maximum while the real part — the in 
phase component — reaches zero. Following, the real 
part changes sign around 1^ [H, HI] . The link with self- 
focusing can then be thought of as follows: the in phase 
component of the space-charge field — the real part — 
is a local photorefractive response, responsible for self- 
focusing or defocusing, whereas the out-of phase one -the 
imaginary part — is responsible for beam bending. This 
phenomenological interpretation has been confirmed by 
an experiment [1^. Later on, analytical derivations and 
numerical simulations [13| came to the same predictions. 
As stated before, however, we were never able to repro- 
duce the experiment, for a reason we will now try to give 
hints of. 

First of all, in the case of InP:Fe, there is a close and 
fairly simple relationship between Id and Ir'- 



Id 



Ir f-lp Cn 



IJ-nCp TlTO 
X . 



PTO 



(9) 



We have evaluated the above ratio to 0.65 with values 
issued from the literature [Slj, [H, HI] , a numerical value 
which is fully compatible with those we have measured 
at 1.56/xm, the discrepancy for the measures at 1.06/im 
being accounted for by the strong absorption impairing 
the Id estimation from figures [5] and [31 

The above ratio ^ evidences the ratio [Fe^+j / [Fe^+] 
in the dark (jito/pto)- Therefore, the oxydo-rcduction 
state plays a significant role, as well as the doping of 
the samples, in the ratio of the two characteristic inten- 
sities. Therefore, a different doping but also a different 
[Fe^+] / [Fe^+] ratio might be the cause of the discrep- 
ancy between our experimental results, which do not ex- 
hibit a crucial role for the resonance intensity, and those 
previously reported. 

However, our theoretical analysis, model II, does not 
either bring the resonance intensity forth. A more univer- 
sal model is in development in order to predict a doping 
dependent behavior and thus unify model II and model 
HI. 



V. CONCLUSION 

In the light of recent experimental results, we have 
shown that despite more complex photorefractive charge 
transport mechanisms, steady-state photorefractive self- 
focusing and solitons in InP:Fe can be described by the 
same model as was done for the simpler one-free-carrier 
model, provided diffusion is neglected. We have evi- 
denced in the process the role of the dark irradiance. Our 
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model does not. however, evidence any crucial role for the 
resonance intensity, a role which is to be the subject of 
future studies. 
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